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Enhanced Spin Susceptibility toward the Charge-Ordering Transition 
in a Two-Dimensional Extended Hubbard Model 

Kazuyoshi Yoshimi*, Takeo Kato and Hideaki Maebashi 

Institute for Solid State Physics, University of Tokyo, Kashiwanoha, Kashiwa, Chiba 277-8581 

Based on the non-skeleton diagrammatic expansion satisfying the compressibility and spin- 
susceptibility sum rules, we investigate static charge and spin responses in a two-dimensional 
extended Hubbard model with the nearest-neighbor Coulomb repulsion in the vicinity of its 
charge-ordering transition point. In this expansion, we can calculate approximate charge and 
spin response functions by systematic inclusion of vertex corrections, from which we obtain 
the uniform susceptibility equal to the so-called q-limit of the response function and the 
second-order transition point as a divergent point in the same response function at some 
finite wave-number vector. It is shown that the reentrant charge-ordering transition, which 
has already been observed in the random-phase approximation (RPA), remains to take place 
even though the vertex corrections are included beyond the RPA. As a prominent effect 
of the vertex corrections, we find that the uniform spin susceptibility is enhanced due to 
charge fluctuations developing toward the charge-ordering transition. We give a qualitative 
comparison of this enhanced spin susceptibility with the experimental results on the quasi- 
two-dimensional organic conductors, together with its explanation in the Landau's Fermi- 
liquid theory. 

KEYWORDS: conserving approximations, vertex corrections, extended Hubbard model, charge 
ordering, organic conductors 



1. Introduction 

Low-dimensional organic conductors exhibit various electronic properties such as super- 
conductivity, magnetism and charge ordering. Among a variety of phase transitions appearing 
in organic conductors, charge-ordering (CO) phenomena have recently attracted much atten- 
tion in the context of strongly correlated electron systems. 1 It has been established that these 
CO phenomena are caused by long-range Coulomb interaction between conduction electrons. 
Quasi-two-dimensional material ^-(BEDT-TTF^RbZm is one of typical compounds most ex- 
tensively studied for CO. This system shows a sharp metal- insulator transition characterized 
by a rapid exponential increase of the resistivity below Tmj = 195K, 2 where a structural tran- 
sition accompanying the doubling of the lattice periodicity along c-axis occurs. Clear evidence 
of charge disproportionation in the CO phase has been reported by NMR and Raman/infrared 
spectrum, whereas long-range ordering has been observed by X-ray diffraction experiment. 3 
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In contrast to remarkable change in the charge degree of freedom, the spin susceptibil- 
ity is frequently less sensitive to the CO transition, but exhibits non-trivial dependence on 
temperature. For example, the spin susceptibility of #-(BEDT-TTF)2RbZn4 has a broad peak 
around the metal-insulator transition temperature Tmi without any sign of singularity. A para- 
magnetic phase is kept down to 30K at which a spin-density-wave (SDW) transition occurs, 
showing the so-called spin-gap behavior. This behavior of the spin susceptibility above 30K is 
analyzed qualitatively by the Bonner-Fisher curve of the Heisenberg model, which is an effec- 
tive model of localized electrons in a strong coupling region. Although this analysis may be 
justified in the insulating phase, it is nontrivial to apply it to the metallic phase above Tmi- It 
is also clear that this effective spin model cannot treat the effect of charge disproportionation 
and its fluctuation properly. Oppositely, if one starts with a metallic state, and incorporates 
the Coulomb interaction as a perturbation, the temperature-dependence of the spin suscepti- 
bility above Tmi is captured by correction to temperature-independent Pauli paramagnetism. 
This correction is expected to be induced by development of quantum fluctuations near the 
CO transition. This weak-coupling approach has not been studied in the previous theoretical 
studies of the CO phenomena. 

So far charge fluctuations developed in the uniform metallic phase near the CO transition 
points have mainly been studied by the random-phase approximation (RPA). In this approxi- 
mation, the CO transition is determined by a divergence in the static charge response function 
at a finite wave-number vector since it indicates instability of the uniform metallic state. If 
the transition is of the second order, the transition point thus determined agrees with the 
one obtained by the self-consistent Hartree approximation, which has succeeded in obtain- 
ing various spatial patterns of CO states in quasi-two-dimensional organic conductors. 4 ' 5 We 
note that the RPA has also been applied to the problem of superconductivity next to the CO 
phase. 6-10 

It is, however, emphasized that the RPA is not suitable for understanding the non-trivial 
temperature-dependence of the uniform spin susceptibility because the spin susceptibility in 
the RPA is not affected at all by enhanced charge fluctuations near the CO transition points. 
In order to include effect of this charge fluctuations, we need to study vertex corrections to 
the response function. As leading vertex corrections, the Maki-Thompson (MT) type and the 
Aslamazov-Larkin (AL) type are well known. In the two-dimensional (2D) Hubbard model 
near the half-filling, these two types of vertex corrections have been investigated, respectively, 
to show that nested antiferromagnetic spin fluctuations lead to the tiny Drude weight 11,12 
and the enhanced charge compressibility. 13 ' 14 

In this paper, we investigate the vertex corrections to the charge and spin response func- 
tions in the 2D extended Hubbard model (EHM) at the 3/4-filling in the presence of the 
nearest-neighbor Coulomb repulsion. We show that the charge fluctuations lead to an en- 
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hancement of the spin susceptibility toward the CO transition at which both of the MT and 
AL type vertex corrections give an important contribution. It is, here, noted that if one takes 
some approximation, one encounters the problem that the compressibility (spin-susceptibility) 
sum rules is not necessarily guaranteed to hold; namely, the uniform charge (spin) suscepti- 
bility defined by the derivative of the electron number (magnetization) with respect to the 
chemical potential (external magnetic field) does not necessarily agree with the (/-limit of the 
charge (spin) response function. To avoid this problem, we construct a theoretical formalism 
based on the non-skeleton diagrammatic conserving approximations, assisted by a set of the 
so-called Hedin's equations, 15 which provide a formal scheme for iterative generation of a 
series of the skeleton diagrammatic conserving approximations. 16 

In our non-skeleton diagrammatic formalism, the compressibility and spin-susceptibility 
sum rules are satisfied at each level of the approximation as in the skeleton diagrammatic con- 
serving approximations, 17 ' 18 although our approximate response function has a rather simple 
form which enables us to make an actual calculation in the lower levels of the approxima- 
tion: The Oth level of the approximation corresponds to the Hartree approximation for the 
self-energy and the RPA for the response function. By iterative generation applied to this 
Oth level, we get the 1st level of the non-skeleton diagrammatic conserving approximation 
abbreviated by 1NSCA, in which the response function includes naturally the MT and AL 
type vertex corrections while the self-energy is described by the RPA. 

This paper is organized as follows. In § 2, we introduce the non-skeleton diagrammatic con- 
serving approximation to formulate systematic inclusion of the vertex corrections. For readers 
who are not interested in details of theoretical formulation, this section may be skipped. In 
§ 3, we apply this formalism to the EHM on the 2D square lattice and then make actual 
calculations of the static charge and spin response functions in the 1NSCA. From the results 
of these response functions, we obtain the spin and charge susceptibilities in the uniform 
metallic phase near the CO transition point. After giving discussions including relation to the 
Fermi-liquid theory and comparison with the experiments in § 4, a summary is given in § 5. 

2. Formulation 

2.1 Exact relations 

In this section, we consider a generic single-band model Hamiltonian H on a lattice with 
band dispersion and electron-electron interaction v af7 i (q) to formulate systematic inclusion 
of vertex corrections: 



fc,CT 




(1) 



a,a' k,k' ,q 
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where Q is the volume of the lattice system, cj^ (cfe CT ) is the creation (annihilation) operator 
of an electron with a wave number vector k and a spin a (= +(— ) for j (I))- Because of 
usefulness for describing the exact relations, we have introduced fi a = ji + oh with /it and 
Ai being a chemical potential and an external magnetic field, respectively. Throughout this 
paper, we take the limit h — ► at the end of calculation. 

The charge and spin densities for the Hamiltonian 7i are given by n = n a and m = 
driu, respectively. Here, n a is the electron number density for a spin a; n a can be related 
to the single-particle Green's function G a (k) by 

n<T = f e u ">G a (k), (2) 

J h 

where J fc denotes (T/Q)^2 k = (T/0,)J2kYli W1 th k being a combined notation of a wave 
number vector k and a fermionic Matsubara frequency iq = (2/ + l)7riT with an integer /, 
and r] is a positive infinitesimal. We can write G a {k) as 

G,(*) = — ~— -, (3) 

iei + Ha- e k - lu a (k) 

where S CT (/c) is the self-energy in which the Hartree contribution is subtracted while the 
chemical potential is shifted by this Hartree contribution as 



A*<r = At. 



a-^v(0)v. (4) 



The charge and spin response functions are related to the density-density response function 
Xaa'(q) as 

WW = g 2 ( 5a ) 



cr,(T' 



where q represents a combined notation of a wave number vector q and a bosonic Matsubara 
frequency \uji = 2W\T with an integer Z; x<tct'(<7) can be written in terms of the one-interaction 
irreducible part Xcc 1 {l) as 

Xaa'(q) = Xaa'(q) ~ ^ Xcro-i (q)v {Q)Xa 2 a'(q), (6a) 

0"1,0"2 

X*o>{q) = ~ [ G a (k + q)G a (k)A aa ,(k;q). (6b) 

J k 

Here, A. ac7 r (k; q) is the vertex function, including the vertex corrections to the response func- 
tions. 

The uniform charge and spin susceptibilities (per spin), which are of central interest in this 
paper, can be defined by Xc = (l/2)(<9n/<9/x) = (1/2) J2acr'(^ n ^/^^^') an d Xs = (l/2)(dm/dh) 
= (V2)E ffff ' (Tcr'firitr/dntr'), respectively. Then Xc and Xs should be equal to the so-called 
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g-limit of the charge and spin response functions as 

1 Oti 

*=ifjr =*«*«»■ < 7b > 

Equations (7a) and (7b) are, respectively, called the compressibility and spin-susceptibility 
sum rules, which hold when both the response functions and the isothermal susceptibilities 
are exactly calculated. If one takes some approximation, however, these equations are not 
necessarily guaranteed to hold. In this paper, we shall take an approximation for the response 
function and determine the charge-ordering transition point such that 

X NN (Q,ty = 00 at some finite q = Q* . (8) 

In order to study the behaviors of the charge and spin susceptibilities near the transition point 
thus determined, it is desirable that eqs. (7) should be satisfied in the approximation we shall 
take. For construction of such an approximation, it is crucial to satisfy the following (/-limit 
Ward identity: 

A CTCT ,(A;;0)=<W-§^ • (9) 

We can show that the sum rules (7) hold automatically in arbitrary approximation satisfying 
the Ward identity (9) as follows: By eqs. (2), (3) and (4), we can easily see 

Substituting eq. (9) into eq. (10) and using eqs. (6), we can, in fact, confirm that eqs. (7) hold. 

Before closing this subsection, we give some comments on a set of the so-called Hedin's 
equations in the quantum many-body theory. 15 In addition to eqs. (6), the following equations 
make up this set: 

t a (k) = -^2 I G a (k + q)V ar7/ (q)A aa ,(k;q), (11a) 
where V aa > (q) represents the effective electron-electron interaction given by 

(71(72 (q)V a2 a'(q)- (Hb) 

"1,(72 

If vertex corrections are neglected, i.e. the vertex function is approximated by A aa i(k; q) = 6 a>a i 
in eqs. (6) and (11), one gets renormalized RPA. It is, however, noted that neither eqs. (7) 
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nor eq. (9) holds in the renormalized RPA. Generic algorithm to include vertex corrections, 
preserving eq. (9) and therefore eqs. (7), was proposed by Baym and Kadanoff (BK) 17 ' 18 
based on the skeleton diagrammatic expansion with respect to the dressed Green's function. 
In principle, the exact self-energy can be obtained from eqs. (6) and (11) combined with this 
BK algorithm. 16 

In the following subsections, we introduce two approximation procedures, both of which 
satisfy the Ward identity (9); one is the well-known procedure based on the mean-field ap- 
proximation in § 2.2 and the other is a non-trivial approximation procedure including vertex 
corrections in § 2.3. For details of a systematic procedure for generating a series of such 
approximations, see Appendix. 



2.2 Mean-field and random-phase approximations 

In the Oth level of the approximation, we neglect S CT (/c) in eq. (3), i.e., f^\k) = 0, leading 
to the Hartree or mean-field approximation for the single-particle Green's function as 

Gf{k) = - 1 . (12) 

iei + Ha-£h 

By virtue of the Ward identity (9), A^ a ,(k;q) should be approximated by 

A^,(fc;<7) = <W- (13) 

Substituting eq. (13) into eq. (6b), we get the one-interaction irreducible response function in 
the Oth level as 

J&{q) = UlK,a' = - f Gf(k + q)G$Xk)5 a , a , . (14) 

Then eq. (6a) leads to the (unrenormalized) RPA for the response functions. This approxi- 
mation is, however, not appropriate to the present study of the uniform susceptibilities; the 
uniform (q = 0) and static (uji = 0) response functions are not affected at all by charge 
fluctuations developed at a finite wave-number vector q = Q* , because the static response 
function is calculated in a closed form for a given q, and has no mode coupling. 



2.3 Inclusion of vertex corrections 

In our approximation procedure, the Ward identity (9) is required to hold as discussed 
in § 2.1. For this purpose, it is essential to utilize the concept of the Hedin's formalism con- 
structed by the skeleton diagrams, 15 but we employ the non-skeleton diagrammatic formalism. 
Substituting eqs. (13) and (14) into eqs. (11), we get the self-energy in the 1st level of the 
approximation as 

tg\k) = -l G< 0) (fc + #i?(g). (15) 

Jq 
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G«(fc) G?\k) 




Fig. 1. The Feynman diagrams (a) for the effective electron-electron interaction V , (q), (b) for the 
approximate single-particle Green's function G<p(k), and (c) for the approximate vertex function 
A^/. The second term of r.h.s. in (c) corresponds to the Maki- Thompson (MT) type vertex 
correction, while the sum of the last two terms corresponds to the Aslamazov-Larkin (AL) type 
vertex corrections (A^, = A A -M + A^/ 2 ). 

V <7<7 <7<J' 1 <7<7 / 



Here V^J (q) is the effective interaction in the RPA defined by 

v£> (q) = v aa , (q) - £ ««i (9)x8U («)^ (?) , (16) 

Cl,(T2 

where the diagram is shown in Fig. 1(a). The approximate Green's function is then given by 

[GW^-^^Wr 1 -^^), (17) 

where the diagram is shown in Fig. 1(b). 

The vertex function A^}\(k;q) carrying a momentum q is constructed by replacing the 
internal G { °)(k') line in (k) by G { °) (k' + q)G < i)(k') in all the possible ways. Here, note 
that all the internal momenta k' are replaced so as to conserve the wave numbers and the 
Matsubara frequencies at all the internal and external vertices. The diagrams for the vertex 
function thus obtained are shown in Fig. 1(c). By this procedure, we can derive various types 
of the Ward identity, one of which is the Ward identity (9) for the g-limit charge vertex. This 
identity can, in fact, be proved by seeing the fact that the differential operation in eq. (9) is 
equivalent to an operation in which one first plucks out the internal G^)(k') line in T,^\k) 
in all the possible ways and then replaces it with G^){p) 2 ; this operation is nothing but the 
vertex insertion described above for q = 0. 

The vertex corrections can be written as the sum of the two contributions as 

A£, (fc; (?) = <V + aMT (^5 ?)<W + Kb (k; q) . (18a) 
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(a) X SC (q) = 



(b)x MT (<?) = 



o 



,~.AL 



(C) Xaa'iQ) = 



€F<o> 



Fig. 2. The Feynman diagrams for the one-interaction irreducible response function (a) without 
vertex corrections, (b) with the Maki-Thompson (MT) type vertex correction and (c) with the 
Aslamazov-Larkin (AL) type vertex corrections. The thin and thick solid lines with arrows rep- 
resent the bare and dressed Green's functions, G^\k) and G^\k), respectively. The thick wavy 
lines represent the effective electron-electron interaction vjfj, (q) . 



The Maki-Thompson (MT) type vertex correction A MT (k;q) is written as 

A MT (k;q)= [ V^(q 1 )G^(k + q 1 + q)G^(k + qi ), (18b) 

while the Aslamazov-Larkin (AL) type vertex correction A^,(k;q) is written as 

A^ L , (k; q)= [ GP (k - qi )V^ ( qi )V^ ( Ql + q) 
J qi 

x [7^(9i; 9) +c.c] , (18c) 

where 7,7(91; q) is the fermion loop with three vertex insertions carrying momenta q\, q and 
-qi-q given by 

la(qi;q) = - [ G^(p)G^(p + q 1 )G^(p + q 1 + q). (18d) 

The one-interaction irreducible part x^J, (q) in the present approximation is given as 

*£,(</) = - / GW(fc + g )GW(fc)A« (k;q). (19) 
There are three contributions in x^,(<z): 

= X S °(q)S^ + X MT (q)S^' + x£(«z) • (20) 

The first one corresponds to a bubble diagram which has only the self-energy correction (SC) 
, and the latter two include vertex corrections. The Feynman diagram for each component is 
shown in Fig. 2. Note that these diagrams are a little different from the usual ones appearing 
in previous literatures; in both x MT (<7) and x^/(g), two °^ tne Fermion lines are the dressed 
Green's functions G^ , but all the other lines are the bare Green's functions G^ ■ In the non- 
skeleton diagrammatic formalism, only this combination guarantees the compressibility and 
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spin-susceptibility sum rules, eqs. (7). In the 1st level of the non-skeleton diagrammatic con- 
serving approximation (1NSCA), the spin and charge response functions are then calculated 
by eqs. (5) and (6) with the approximate one- interaction irreducible function x^L(q). These 
response functions contain the leading contribution from nontrivial vertex corrections and 
introduce mode coupling between charge and spin fluctuations with different wave numbers. 
Since the MT and AL type vertex corrections have a crucial role near a second-order transition 
point, we can obtain physically comprehensive results at this level of approximation. In the 
next section, we study a role of these vertex corrections near the CO transition based on this 
approximation . 

Finally, we mention an extension of the present approximation. We can construct the next 
level of approximation by making the next self-energy from the lastly obtained vertex and 
response functions. Performing this procedure iteratively, exact series of diagrams for all the 
functions are generated as in the skeleton-diagrammatic expansion, 16 although it is practically 
difficult to sum up them. Details of this formal extension are given in Appendix. 

3. Numerical Results and Discussions 

3.1 Model 

In this section, we consider the single-band (extended) Hubbard model on the 2D square 
lattice. In order to investigate effect of fluctuations on the response functions near the CO 
transition point, we include the nearest-neighbor Coulomb interaction V into the model in 
addition to the on-site Hubbard interaction U. The Hamiltonian of the single-band EHM is 
given by 

fc,<7 1,(7 (ij) 

where the band dispersion is assumed to be = 2t[cos(k x a) + cos(k y a)} with t and a being 
the transfer integral between nearest-neighbor sites and the lattice constant, respectively. The 
chemical potential fi is determined to fix n = ^2 a n a = 3/2 by using Eq. (2) and a denotes 
the spin with a direction opposite to a. In the following, we set t = 1, a = 1 and &b = 1. 

The Hamiltonian Hehm corresponds to a special case of the generic single-band model 
Hamiltonian (1) with the bare Coulomb interactions given by 

vM = V(q), (22a) 

v a a(q) = U + V(q), (22b) 

where V(q) = 2y(cosg :r + cos q y ). For the charge and spin channels, the bare interactions are 
represented by v c (q) = v aa (q) + v a9 {q) = U + 2V(q) and v s (q) = v aa (q) - v a& {q) = -U, 
respectively. Then the charge and spin response functions can be written in the form 
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Fig. 3. Phase diagram on the (V, T) plane at U — 3. The dashed and solid curves represent the CO 
transition curves in the RPA and in the 1NSCA, respectively. 



1 + M9)Xs,s,(?) 



Where X NN (q) = Xaa(q) + Xaa(l) and X Sz S z (<?) = Xaa(q) - Xaail)- 

In the RPA described in § 2.2, x NN (q) and x SzSz (q) are approximated by Xy N (q) = 
Xs s (^) = Xo{q)i where Xo(q) 1S defined by eq. (14). In the 1st level of the non-skeleton 
diagrammatic conserving approximation (1NSCA) described in § 2.3, they are approximated 
by 

xSUq) = X SC (q) + X MT (q) + X^(q) + x£(<f), (24a) 

41 («) = * S °(9) + X^fo) + X^fo) " ( 24b ) 
where X SC (<?) is a contribution associated with a bubble diagram (Fig. 2(a)) which has only 
the self-energy correction (SC) while x MT {q) an d xwt'i^) i ncm de the Maki-Thompson (MT) 
and Aslamazov-Larkin (AL) type vertex corrections (Figs. 2(b) and 2(c)), respectively. Note 
that the difference between the one-interaction irreducible response functions x^ N (q) an d 
Xg 1 5 (q) comes only from the AL contribution. 

In order to calculate the convolution form, we use the Fast-Fourier- Transform (FFT) 
algorithm. The first Brillouin zone is divided into 64 x 64 meshes. The frequency sum is 
terminated at e c whose value is about 40 times as large as the band-width W = 8t for 
T = 0.1. In this section, we first determine the phase diagram on the (V, T) plane. Next, we 
calculate the ^-dependence for the uniform spin and charge susceptibilities. In the following, 
we fix U = 3, at which a SDW instability is not observed for T > 0.05. 
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Fig. 4. Static charge response function in the 1NSCA for (T, V) = (0.05, 1.7), (0.1, 1.7), (0.2, 1.6) 
and (0.3, 1.5) at U = 3. The peak position of the response function continuously shifts from 
Qcb = (7r,7r) to Q ic « (37r/4,37r/4) for T < 0.3. 

3.2 Phase diagram 

The phase diagram on the (V, T) plane at £/ = 3 is shown in Fig. 3, where the CO 
transition curves are determined by the divergence of the static charge response function, i.e., 
by l + v c (q)x NN (q, 0) = in the RPA (the dashed curve) and in the 1NSCA (the solid curve). 
Since the region of the CO state in the 1NSCA is smaller than that in the RPA, we can 
understand that the CO transition is suppressed by the charge fluctuations. It is, however, 
noted that the overall shapes of the transition curves in these approximations are similar to 
each other. In particular, the reentrant transition (metallic — > CO — > metallic with decreasing 
temperature) are observed both in the RPA 6, 19 and in the 1NSCA; namely, the reentrant CO 
behavior remains even if the vertex corrections are taken into account. 20 

In Fig. 4, we show the static charge response function in the 1NSCA for the uniform 
metallic state near the CO transition curve in the phase diagram. At high temperatures, the 
q-dependence of x NN (?> 0) does not depend on q very much, resulting in checkerboard-type 
charge ordering at q = Q c b = (ft, 7r) at which v c (q) takes a minimum value. With decreasing 
temperature, the existence of the Fermi surface causes the strong q-dependence of x NN (q), 
so that the peak position of the static charge response function changes from q = to 
q = Qic ~ (3vr/4, 37r/4) for T < 0.2, leading to an incommensurate CO instability at q = Q lc . 
The reentrant transition observed in Fig. 3 is due to nonmonotonic temperature-dependence 
of Xjvjv(Qcb)O) with a fixed V; namely, as temperature increses, X NN {Qcb,fy once increases 
until T < 0.5 but decreases at high temperatures. This unusual increase of Xivjv(Qcb)O) for 
low temperatures originates from broadening of the incommensurate peak of Xjvjv(9'0) at 

q = Qic 
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3.3 Uniform susceptibilities 

Let us move on the main result of this paper. In our approximations, the compressibility 
and spin-susceptibility sum rules (7) are satisfied, so that the uniform charge and spin suscep- 
tibilities are equal to the g-limit of the response functions as Xc = Xnn(®) an d Xs = Xs z s z (0)> 
respectively. In Fig. 5, we depict the results for Xs/xo (thick solid curve) and Xc/xo (thin 
solid curve) in the 1NSCA for T = 0.1 and U = 3 where the non-interacting susceptibility 
is denoted by xo = Xo(0). In the same figure, we also represent the results for Xs/xo (thick 
dashed line) and Xc/xo (thin dashed curve) in the RPA for comparison. Note that the CO 
transition occurs at V = 1.788 in the 1NSCA for T = 0.1 and U = 3. 

From Fig. 5, we find that the spin susceptibility Xs calculated by the 1NSCA is enhanced 
toward the CO transition point of V = 1.788. Since Xs is independent of V in the RPA, we 
can conclude that this enhancement is caused by the vertex corrections. On the other hand, 
the charge susceptibility Xc decreases with the increase of V both in the RPA and in the 
1NSCA; the effect of vertex corrections on Xc is small as indicated by the fact that there is 
only a slight difference between the results of Xc in these two approximations. 

For detailed explanations on these ^-dependences of Xc an d Xs, we rewrite the g-limit 
of the self-energy, MT and AL contributions in eq (24) as Xc C = Xs = X. (0), xf T = 

~MT _ ~MT (0); ~AL = ~AL (0) + ~AL (Q) and ~AL = ~AL (Q) _ ~A_L (0) for simplicity . In the 

following subsections, we will investigate each of these contributions to the spin and charge 
susceptibilities individually. 




Fig. 5. Uniform spin and charge susceptibilities x s and Xc divided by the non- interacting suscepti- 
bility xo versus the nearest-neighbor Coulomb repulsion V for T — 0.1 and U — 3. Xs/xo an d 
Xc/xo are plotted by thick and thin solid curves, respectively, in the 1NSCA and by thick and 
thin dashed curves in the RPA. Note that the incommensurate CO transition occurs at V = 1.788 
in the 1NSCA. 
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Fig. 6. The ^-dependences of xf C [= Xc°] 



xf h and x^ L shifted by their values at 



V = 0. The temperature and the on-site Coulomb interaction are fixed as T = 0.1 and U = 3, 
respectively, where the CO transition occurs at V = 1.788. 



3.3.1 Spin susceptibility 

The uniform spin susceptibility is calculated in the 1NSCA as 

Xs = -jr. ; , (25) 

where X s s (0) = X^ C + X^ T + X^ L - From eq. (25), it is clear that the Independence of Xs 
comes only from x^\ (0)- Therefore, for a positive U, the enhancement of Xs observed in 
Fig. 5 is caused by the increase of Xy s (0). In Fig. 6, the Independence of Xs C , X^ T and xf L 
are shown by subtracting their values at V = 0, respectively. It is found that, as V increases, 
X^ L and x^ T increase, while xf C decreases. 

First, we discuss the behavior of Xs and X^ T - Up to the leading correction with respect 
to U and V ', they can be written separately as 

X S s C = Xo + Xu C + Xv C > (26a) 
xf T = x¥ T + xT- (26b) 

Since we are interested in the ^-dependence of the spin susceptibility, we only focus on X\F 
and Xy T > which are proportional to V. At low temperatures, these corrections are expressed 
by the integral of the interaction on the Fermi surface Sf as 21 

X S y° = ~XT (27a) 

xT = f dkdk'V(k - k')5(e k - p)5{E k , - M ) 

= (V(q))s F . (27b) 
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Fig. 7. The n-dependence of the average of the nearest-neighbor Coulomb interaction taken over the 
Fermi surface: The inset shows the n-dependence of the Fermi surface for the non-interacting 
system. 



The sign of (V(q)) depends on both V{q) and the shape of the Fermi surface. In the present 
simple model with nearest-neighbor Coulomb interaction on the square lattice, the potential 
V(q) = 2V{cosq x + cosq y ) is positive for small momentum transfer around q ~ 0, and 
negative for large momentum transfer around q ~ (ir, it). We show {V(q))s F for 1 < n < 2 in 
Fig. 7, together with the change of the Fermi surface in the inset. At the half filling (n = 1), 
(V(q)) becomes zero because positive contribution of small momentum transfer is canceled by 
negative one of large momentum transfer. As the filling increases, the Fermi surface becomes 
smaller, and therefore negative component is reduced. As a result, (V(q)) becomes positive, 
and increases until n ~ 1.9. Finally, (V(q)) rapidly decreases toward n = 2 where the Fermi 
surface disappears. We note that the above discussion is also applicable for < n < 1 by 
replacing n with 2 — n. Thus, at least in the leading correction with respect to V, the MT 
contribution Xy T ^ s positive for the 3/4-filling (n = 3/2), leading to the enhancement in the 
spin susceptibility. It is, however, noted that by eq. (27a), this enhancement due to the MT 
contribution is canceled by the self-energy contribution X\F for small V. These features of 
X\F and X\/ T are consistent with the results shown in Fig. 6. 

Next, we discuss m terms of the effective interaction V^ A (q) = V^}(q,0). Since 
the MT contribution is almost canceled by the self-energy one even near the CO transition 
point of V = 1.788 as seen in Fig. 6, the net enhancement of the spin susceptibility toward 
the CO transition results mainly from x^ L - I n the AL type diagrams shown in Fig. 2(c), 
the two of the effective interactions (the thick wavy lines) carry the same momenta for the 
uniform {q = 0) susceptibility and their static components around the ordering vector Q* are 
dominant near the CO transition point; on the other hand, the Green's functions themselves 
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(the thick or thin solid lines) are less sensitive to the increase of charge fluctuations and give 
positive contribution to x AL . We can then make a rough estimation of x AL as 

x AL oci5>V^ A (Q*) 2 



<J,(J' 



= v c KPA (Q*)v s KPA (Ql, (28) 

where V c RPA (g) and V RPA (q) are given by 

V c RPA (q) = - V ;[ Q ] . - , (29a) 

l + v c (q) X0 (q,0y 

V s RPA (q) = V ; {q \ - . (29b) 

It is noted that in the paramagnetic phase without charge ordering, the denominators in 
eqs. (29), i.e., 1 + v c (Q*)xo(Q*, 0) and 1 + v s (Q*)xo(Q* , 0) are both positive. Because charge 
fluctuations are developed only for the wave-numbers satisfying v c (Q*) = U + 2V(Q*) < 0, we 
can then see that amplitude of V RPA (Q*) is developed toward the CO transition point keeping 
its sign negative. On the other hand, V RPA (Q*) is always negative because v s (Q*) = — U < 
for a positive U. Therefore, eq. (28) indicates that x AL increases toward the CO transition; 
this behavior is consistent with the result in Fig. 6. 

In conclusion, charge fluctuations developed near the CO transition point make the uni- 
form spin susceptibility increase through the vertex corrections; the behavior of Xs in Fig- 5 
can, in fact, be understood qualitatively by the increase of x 



-AL 



3.3.2 Charge susceptibility 



The uniform charge susceptibility is calculated in the 1NSCA as 

Xc = (1) ; - / , (30) 

where X^(0) = Xc° + X^ T + X. AL - The overall ^-dependence of Xc is dominated by 8V, 
i.e., the last term in the denominator in eq. (30), so that Xc decreases against the increase 
of V as we have seen in Fig. 5. It is, however, interesting to consider the effect of charge 
fluctuations on the irreducible charge susceptibility X^(0). This consideration may suggest 
that Xc increases eventually towards the CO transition if the higher-order vertex corrections 
are properly taken into account beyond the 1NSCA. 

Noting that Xc° = Xs° an d = xY' 1 \ we see that the self-energy and MT contributions 
almost cancel each other in the charge susceptibility as well as in the spin susceptibility dis- 
cussed in the previous subsection; an important difference between the irreducible charge and 
spin susceptibilities is then ascribed to the AL contribution. By using the same approximation 
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in the previous subsection, we can make a rough estimation of x AL as 

= [V C RPA (Q1] 2 + [V S RPA (Q*)} 2 . (31) 

In contrast to x AL dependent linearly on F C RPA (Q*), eq. (31) indicates that x AL increases 
toward the CO transition in proportion to the square of V^ PA (Q*). Therefore, Xc is expected 
to increase when x AL changes considerably to overcome the increase of 8V in eq. (30). In 
order to study this possibility in details, charge fluctuations should be fully treated beyond 
the present approximation. This problem will be studied elsewhere. 22 

4. Discussions 



4-1 Relation to the Fermi-liquid theory 

In this paper, the spin and charge susceptibilities are studied up to the leading vertex cor- 
rections with respect to charge fluctuations. We expect that the present result gives a correct 
behavior of the susceptibilities toward the CO transition until the charge fluctuations are not 
developed so much. Just near the transition, however, the higher-order vertex corrections may 
become relevant due to strong charge fluctuations. Although the treatment of the higher-order 
corrections is beyond the scope of this paper, it is feasible to summarize its effect in terms of 
the Fermi-liquid theory, which is a natural description of the metallic phase. 

Following a standard microscopic derivation of the isotropic Fermi liquid, 23 the uniform 
susceptibilities may be written as 

x NN (Q) _ p (32) 

Xc ~ l + v c (0)x NN (0) "1 + P/ S ' ( j 

Y - Xs z s z (°) _ P to 2h \ 

Xs l + v s (0)x SzSz (0) 1 + [62b) 

where /q and fg are the charge and spin channels of the Landau's quasi-particle interaction, 

respectively, and p is the renormalized density of states at the Fermi level in proportion to 

the effective mass of the quasi-particle. For comparison with the results in the last section, 

it is useful to introduce the "screened" quasi-particle interaction in the mean field by /q = 

/* - Uc (0) = fg - U - 8V and f§ = / a - v s (0) = fg + U. We can then write the irreducible 

charge and spin susceptibilities as 

*»» (0) = TT^ (33a) 

WO) = 03b) 
For |p/q| <C 1 and \pfo \ <C 1, in particular, we can expand these irreducible susceptibilities as 
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f MT 
JO 


fs,AL 
JO 


?a,AL 
J 


(a) charge: v c (Q*) < 








(b) spin: v c (Q*) > 


+ 




+ 



Table I. Trends toward an increase (+) and a decrease (— ) in the "screened" quasi-particle interaction 
(a) for enhanced charge fluctuations, and (b) for enhanced spin fluctuations, based on the single- 
band EHM. Note that v c (Q*) is necessarily negative in (a), but v c (Q*) is assumed to be positive 
in (b), with Q* being the ordering vector at which either charge or spin fluctuations are developed. 

follows: 

X NN (0) « p(l - PU) = p[l - P(fr + /o' AL )L (34a) 
X SzSz (0) » p(l - p/ a ) = p[l - p(/ MT + /^ AL )]. (34b) 

In the above expansion, noting that the self-energy correction is included in p, we have de- 
composed /o and /o into contributions from the MT and AL type vertex corrections as 

fS = /o MT + /o' AL and fS = /o MT + /o* AL , respectively. 

In Table 1(a), we summarize the effects of the leading vertex corrections on the uniform 
susceptibilities in terms of f^, /q AL and /q' AL by comparing the results in § 3.3 to eqs. (34) 
for enhanced charge fluctuations. In Table 1(b), in order to emphasize a characteristic change 
due to these charge fluctuations, we also show a change in the same quantities when not 
charge but spin fluctuations are developed near a SDW transition point. In the latter case, 
the Fermi- liquid corrections have already been discussed in details, based on the Hubbard 
model. 24 Note that in the present calculation for the EHM, the spin and charge fluctuations 
can be controlled simply by changing U and V, respectively. 

From Table I, we see that both charge and spin fluctuations make /q' AL decrease in 
eq. (34a) . This is clear from the fact that x AL increases with development of either charge or 
spin fluctuations as in eq. (31). The important point is that and /g' AL in eq. (34b) tend 
to decrease for enhanced charge fluctuations, but to increase for enhanced spin fluctuations; 
in other words, the charge and spin fluctuations provide attractive and repulsive parts in the 
spin channel of the quasi-particle interaction, respectively. The decrease in _/ u vIT results from 
the exchange effect due to the nearest-neighbor Coulomb interaction V, which is absent in 
the Hubbard model, as we have seen in eq. (27b). But this decrease tends to be canceled 
by the decrease in p originating from the self-energy correction, so that the increase of \s is 
controlled by Jq' AL . 

For /o' AL > ec l- (28) is available not only when charge fluctuations are developed but also 
when spin fluctuations are developed. Noting that charge (spin) channel of the effective inter- 
action is attractive near the CO (SDW) transition point, it is not hard to see from this equation 
that the sign of /q ' AL coincides with that of the spin (charge) channel of the interaction. For 
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enhanced charge fluctuations (V^. RPA (Q*) < 0), K, RPA (Q*) is necessarily negative because of 
v s (Q*) = — U < 0, so that f$ decreases. For enhanced spin fluctuations (V^ PA (Q*) < 0), 
on the other hand, V 6 RPA (Q*) is positive unless v c (Q*) = U + 2V(Q*) < 0, so that /g' AL 
increases. This observation of /q ,AL based on the EHM leads us to a general trend that charge 
fluctuations make Xs increase, while spin fluctuations make Xs decrease. 

4-2 Comparison with experiments 

The result obtained in this paper may be compared qualitatively with experiments on 
the spin susceptibility of quasi- two-dimensional organic materials showing CO. Most relevant 
materials are (9-(BEDT-TTF) 2 RbZn 4 and /3-(meso-DMBEDT-TTF) 6 PF 6 . As mentioned in 
§ 1, the spin susceptibility of the former material increases as temperature is lowered from the 
room temperature to the metal-insulator transition temperature Tmi = 195K; its value changes 
gradually from X s(T = 300K) = 6 x l(T 4 emu moP 1 to X s(T ~ Tmi) = 7 x l(T 4 emu mol -1 . 2 
Similarly, the spin susceptibility of the latter material increases from Xs(T = 300K) = 
6 x lCT 4 emu mor 1 to Xs(T = T M i = 90K) = 1.1 x lCT 3 emu mol" 1 with decreasing tem- 
perature. 31 These behaviors in the metallic phase are consistent with our theoretical result. 
For more quantitative comparison, however, we need to consider several effects neglected in 
the present calculation such as electron-phonon interaction, inclusion of realistic band disper- 
sion, and long-range Coulomb interaction beyond nearest-neighbor molecules. In particular, 
we have to note that the present calculation is based on the weak-coupling approach start- 
ing with a noninteracting metallic state; strong electron correlation, which cannot be treated 
in the present approach, may affect the behavior of the spin susceptibility. Nevertheless, the 
present calculation gives one of possible explanations to enhancement of the spin susceptibility 
in terms of charge fluctuations developed toward the CO transition. 

It is comprehensive to consider the spin susceptibility of k-(ET) 2 X and /3-(ET)2l3, which 
are composed of two-dimensional conducting ET molecules forming dimers. 25 These mate- 
rials may be described by the Hubbard model at the half-filling (one hole per dimer) due 
to molecular dimerization. Therefore, they show a metal-insulator transition by band-width 
control with applying pressure or replacing anions. Antiferromagnetic order appears for low 
temperatures in the insulating side, while superconductivity shows up in the metallic side. A 
most metallic material /3-(ET) 2 l~3 shows a nearly temperature-independent spin susceptibility 
( X s ~ 4 x 10~ 4 emu mol -1 ), i.e., the Pauli paramagnetism from the room temperature down to 
the superconducting transition at T = 3.6K. k-(ET) 2 Cu(NCS) 2 and K-(ET) 2 Cu[N(CN) 2 ]Br 
with a narrower band width are closer to the antiferromagnetic phase in the phase diagram, 
but still in the metallic side. The spin susceptibility of these two materials takes a similar value 
as /3-(ET) 2 i3 for high temperatures, but it decreases as temperature is lowered. This exper- 
imental observation is consistent with our expectation that the increase of spin fluctuations 
make the spin susceptibility decrease as discussed in the previous subsection. 
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One advantage of our approach is that two opposite effects on the spin susceptibility by 
development of charge and spin fluctuations can be understood naturally in the weak-coupling 
approach. We conjecture that this result is valid for many organic materials with a simple 
band structure unless higher-order vertex corrections become relevant just near the transition 
point. 

Finally, we shall give brief comments on other CO materials. (DI-DCNQI)2-M (M=Li and 
Ag) and (TMTTF^A" (X =monovalent anion) are famous quasi-one-dimensional conductors 
showing CO, 1,3 but are out of the scope of the present study because low-energy excitations 
in these systems may be described by the Tomonaga-Luttinger liquid in which spin and 
charge degrees of freedom are separated. a-(BEDT-TTF)2i3 is also a well-known CO organic 
conductor with a quasi-two-dimensional crystal structure. 3 The spin susceptibility of this 
material decreases as temperature is lowered toward the CO transition point at T = 135K. 26 
We conjecture that this behavior apparently opposite to our result is related to its multi- 
band structure with four small Fermi surfaces. 2 The present theoretical discussions on the 
single-band EHM can be generalized straightforwardly to a multi-band model by replacing the 
interaction with an interaction matrix. This extension is, however, non-trivial to be solved in 
a future problem. The CO phenomena have also been observed in several inorganic materials. 
For example, the vanadium bronze /3-Nao.33V20s shows the metal-insulator transition at T = 
135 K accompanied by CO. 28-30 The spin susceptibility of this material is enhanced toward 
the metal-insulator transition with decreasing temperature as indicated by our theoretical 
result for a simple one-band model, although our weak-coupling approach may be insufficient 
for this material with relatively large Coulomb interaction. 

5. Summary 

In this paper, we developed the non-skeleton diagrammatic expansion to satisfy the com- 
pressibility and spin-susceptibility sum rules. Based on this expansion, the leading vertex 
corrections to the static response functions have been examined for the two-dimensional ex- 
tended Hubbard model in the vicinity of its charge-ordering transition point. It was shown 
that the charge-ordering transition is suppressed by charge fluctuations with the overall shape 
of its transition curve almost unchanged. We found that the same charge fluctuations make 
the uniform spin susceptibility increase due to interesting coupling of the charge and spin 
degrees of freedom through the vertex corrections. 

The increase in the spin susceptibility originates from the Maki-Tompson and Aslamazov- 
Larkin type vertex corrections. It is, however, noted that the former contribution is almost 
canceled by the self-energy correction, so that the net enhancement of the spin susceptibility 
can be controlled by the latter contribution. A physical interpretation of this result was also 
given in the Fermi-liquid theory: When charge fluctuations develop toward the transition, the 
effective mass of the quasi-particle tends to decrease but an attractive part in the spin channel 
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of the quasi-particle interaction tends to increase. Since the latter effect overcomes the former, 
the uniform spin susceptibility gets enhanced. 

Our result of the enhanced spin susceptibility towards the charge-ordering transition is con- 
sistent with experiments on 6>-(BEDT-TTF) 2 RbZn 4 2 and /3-(meso-DMBEDT-TTF) 2 PF 6 . 31 
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Appendix: Systematic inclusion of vertex corrections 

In this Appendix, we present an iterative approximation scheme for systematic inclusion 
of vertex corrections to the charge and spin response functions. We construct (A) the single- 
particle Green's function G&\k), (B) the vertex function A.^\(k;q), and (C) the response 
functions x$ n {q) an d X§\ (?)■ These functions are constructed so that the Ward identity 
holds at each level of approximation assigned by an integer i, and the next level of approxima- 
tion assigned by i + 1 is constructed from the ith approximation. Our algorithm of generating 
a series of approximations is almost the same as the one proposed by Takada. 16 Our scheme 
differs from Takada's one on three points; (1) we start with the mean-field Green's function 
as the initial (i = 0) approximation, (2) we use the non-skeleton diagrammatic analysis for 
fla\k) with respect to the combination of dressed Green's function, *\k) and G&\k), 
and (3) we simplify iterative procedure by avoiding direct integral for the Bethe-Salpeter 
equations. Our scheme is suitable to study the leading contribution of the vertex corrections. 

(A) The single-particle Green's function: Suppose the self-energy S^(fe) in which the 
Hartree contribution is subtracted is given as a functional of the mean-field Green's functions 
(a = ±), i.e. 

^ (fc ) = EW (fc:[G £>]). (Al) 

We take the initial function as E^(fc) = (the Hartree approximation). Let the single-particle 
Green's function be 

G«(*;) = [G( r )(fc)- 1 -EW(fc)]-i, (A . 2) 
where the corresponding diagram is shown in Fig. A-l (a). G&\k) is a function of ji^; the 
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Fig. A-l. The Feynman diagrams of the present iterative approximation: (a) the Dyson equation, (b) 
construction of the vertex function, (c) the effective interaction, and (d) the next definition of the 
self-energy. 



particle number and chemical potential are calculated as 

n$= Jj^G${k), (A3) 

/,« = /iW + ^w(0)n«. (A-4) 

a' 

In this paper, we fix fiy = i^a during the calculation, and determine fl a by n£ =3/4 from a 
jla-n curve. 

(B) The vertex function: Following the Kadanoff-Baym scheme, the vertex function 
A^ a ,(k;q) = A^ a ,(k;q : [G± ]) is given as the sum of 5 aa ' and all the possible diagrams 
obtained by inserting the external vertex carrying a momentum q into an arbitrary G^-line 
with a spin a' in each diagram for l (k : [G± ]). This procedure is schematically shown 
by the Feynman diagrams in Fig. A-l (b). Note that the wave numbers and the Matsubara 
frequencies at all the internal and external vertexes are conserved. The vertex function thus 
obtained satisfies various types of the Ward identity, one of which is given by 

A«(*;0)=O-^^. (A-5) 

This identity is proved by the same way as in § 2.3 as follows: Since Y>a\k) depends on fi± 
only through G±, the differential operation on (k) with respect to jj, a corresponds to 
an operation to pick up one internal G^) (p) line in T,a\k) in all the possible ways, and to 
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replaces it with G^)(p) 2 . This differential operation is nothing but the vertex insertion for 
q = as seen in Fig. A-l (b). 

(C) The response functions: The response function xfl'il) * s written in terms of the one- 
interaction irreducible part X^,(q) at the ith level of approximation as 

xt (q) = X®, (q) ~ £ {q)v™ (<?)x^ (q) , (A-6a) 



01,(72 



xt(q) = ~ [ Gf{k + q )G${k)kf a ,{k- q ). (A-6b) 
Jk 

Then, the charge and spin response functions x$ N (q) anci X$\ (q) are calculated as 



A s z s z v 2 



(s) = ^E™'^(<?)- ( A - 7b ) 



As proved in § 2.1, the Ward identity (A-5) leads to the compressibility and spin-susceptibility 
sum rules at the ith level of approximation as follows: 

^^jr-xS-W. (A ' 8a) 
xS" = ^ = xg fc (o)- (a*) 

(D) The self-energy in the next level of approximation: For a systematic improvement in 
the approximation, we require that the vertex functions should be taken into account in a 
manner consistent with the single-particle Green's function. This requirement makes us to 
select the self-energy in the next level of approximation as 

ftj+^k) = "E / G®(k + q)V®{q)A®{k ]q ), (A-9a) 

where the effective interaction vj l J, (q) is given by 

V®, (q) = v aa , (q) - £ v aai (q)x« CT2 {q)V^ (q) . (A-9b) 

The diagrams of these equations are shown in Fig. A-l (c) and (d). With Y,a +1 \k) = f^ +1 \k : 
\G^±\), we can construct the approximation for X^^Hq) an d Xg^Ho) following the processes 
in (A), (B) and (C). 

In the text, we have used two approximation procedures described in § 2.2 and in § 2.3. In 
the present iterative approximation scheme, the former corresponds to i = 0, and the latter to 
i = 1. In principle, we can continue the iterative process as one hopes, although it is difficult 
to continue it for i > 2 practically. 

We note that the approximate response functions approach the exact ones in the limit of 
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i — > oo as 

)hnx®„(q)=X„M, ( A -10a) 

lim X W {q)= Xa zS M. (A-10b) 
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